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Abstract 

The existence of cycles of the matrix maps in Fibonacci class of lattices is 
well established. We show that such cycles are intimately connected with the 
presence of interesting positional correlations among the constituent 'atoms' in 
a one dimensional quasiperiodic lattice. We particularly address the transfer 
model of the classic golden mean Fibonacci chain where a six cycle of the 
full matrix map exists at the centre of the spectrum [Kohmoto et al, Phys. 
Rev. B 35, 1020 (1987)], and for which no simple physical picture has so 
far been provided, to the best of our knowledge. In addition, we show that 
our prescription leads to a determination of other energy values for a mixed 
model of the Fibonacci chain, for which the full matrix map may have similar 
cyclic behaviour. Apart from the standard transfer-model of a golden mean 
Fibonacci chain, we address a variant of it and the silver mean lattice, where 
the existence of four cycles of the matrix map is already known to exist. The 
underlying positional correlations for all such cases are discussed in details. 
PACS No.s: 71. 23. An and 71.23.Ft 
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One dimensional quasiper iodic lattices have been the subject of extensive research in the 
past several years, and we have already gathered a wealth of knowledge about these sys- 
tems [1-11,13-16, and references therein]. Theoretical studies of the spectral properties of 
one dimensional quasiperiodic lattices in particular, have occupied a large portion of all 
these works [see 1-11, for example]. Perhaps the most widely studied model is based on 
the Fibonacci sequence which started drawing interest after the publication of a couple of 
pioneering papers by Kohmoto et al [1] and Ostlund et al [2] . The spectral properties of one 
dimensional Fibonacci chain are exotic. The single particle eigenstates are neither extended, 
nor localized in a conventional sense, and the spectrum is a Cantor set with zero Lebesgue 
measure [1-3]. The existence of self similarity and scaling behaviour of the wave functions 
are two other major properties which mark the difference of such a quasicrystalline chain 
from a periodic or a randomly disordered lattice. A Fibonacci chain consists of two letters 
A and B. The entire sequence is generated by successive application of the substitution 
rule, A -> AB and B -> A. The first few generations are, G = A, G 1 = AB, G 3 = ABA, 
G 4 = ABAAB, and so on, where the letter G indicates 'generation'. The letters A and B 
may denote two different 'atoms' (site model) or, two different 'bonds' separating identical 
atoms (transfer model) [7]. The electronic properties of such systems are well described 
within the standard tight binding formalism, where the Hamiltonian, in the Wannier basis, 
is given by, 

^ = E c *l < ><i|+E[*«l < >0'l+*iiU>(i|] (!) 

i ij 

Here, is the site energy corresponding to the i-th site and Uj = tji is the hopping integral 
between the i—th and the j— the sites. To solve for the eigenfunctions it is useful to find the 
solutions of a set of difference equations which, in the nearest neighbor approximation look 
like, 

(E - €i)ipi = ti ti+ ii) i+ i + £ M _iV>j-i (2) 

Here, ipi is the amplitude of the wave function on the i—th site. The above equation can 
equivalently be cast into a matrix form: 
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which relates the amplitudes of the wave function at three consecutive sites on the chain, is 
called the transfer matrix. 

For investigating the spectral properties of a quasiperiodic chain, one major approach, using 
a method of the dynamical systems analysis, has been to generate first a finite segment (a 
'rational approximant' [1,2,7]) of the Fibonacci chain following the substitution rule. The 
infinite chain is then constructed by periodically repeating this segment. In this way one can 
use higher and higher order rational approximants to attain a true quasiperiodic Fibonacci 
chain in the thermodynamic limit. For any generation I one can construct a transfer matrix 
Mi which is the product of the individual transfer matrices across each atomic site in that 
l—th generation, taken in proper order, i.e., Mj = YiiMi. This product matrix Mi plays a 
vital role in determining the character of the spectrum and the eigenfunctions. The allowed 
eigenvalues for the infinite sequence are extracted from a study of the evolution of the trace 
of the above transfer matrix at successive generations of the Fibonacci sequence [1,2,7]. Two 
important properties of this product matrix for any generation are now well established. We 
refer to them with respect to the golden mean Fibonacci chain, which will play a central 
part in this communication. 

(i) The Trace Map: Following the prescription that generates the sequence, it is easy to see 
that the total transfer matrix for any / + 1-th generation sequence is related to the matrices 
corresponding to two previous generations by the relation, M t+1 = M^Mi. This leads to a 
recursion relation between the traces of M n , viz, 



Xi +1 = IXiXi-x - Xi-2 



(4) 



Here xi = TV(M;)/2. The "allowed" energy values are those for which \xi\ < 1. As the 
generation increases we find that almost all the energy values belong to the escaping orbits 
of the map [1,2], which speaks for a true singular continuous character of the spectrum. 

(ii) The Matrix Map: In addition to the above property, the product matrix Mi at any 
generation / is found to exhibit a cyclic behavior [1-7], viz, Mi +m = Mi for certain energy 
eigenvalues that depend on the model considered. To clarify, let us give the example of the 
so called 'transfer model' of the standard golden mean Fibonacci chain, where the letters A 
and B stand for two bonds, 'Long' (L) and 'Short' (S) respectively. All the on site potentials 
are assumed to be same, and put equal to zero. In this model one finds that, at E — which 
happens to be the centre of the spectrum in this case, the matrix product M; has a six 
cycle, viz, M/ +6 = Mi for I > 2. It may be noticed that such cyclic behaviour of the transfer 
matrix is also common to quasiperiodic chains other than the golden mean [10]. It is this 
cyclic property of the matrix Mj, that has stimulated the present work. The motivation will 
be clear from the discussion that follows. 

The cyclic behaviour of the transfer matrices, as it appears from the previous works in 
general, has been shown to exist for various quasiperiodic chains without any underlying 
physical reason. It is not straightforward to prove the existence of cycles of the matrix 
map for any arbitrary substitutional^ generated sequence. Not only that, even if it exists, 
it is extremely difficult to locate the precise energy eigenvalue for a lattice such as the 
Fibonacci chain in which the spectrum is a Cantor set having a gap at the vicinity of every 
energy, at which such a cycle will take place. Precisely, this is the point that motivates 
us to undertake the present study. We wish to examine whether it is at all possible to 
provide a simple physical picture that will lead to an understanding the origin of cycles of 
the matrix map based on the geometric structure of the quasiperiodic chain concerned. In 
this communication we begin by examining the well known six cycle behaviour at E — 
in the simplest transfer model of a golden mean Fibonacci chain. We show that there 
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exists hidden dimer- like correlations between the constituent 'atoms'. The dimers are not 
explicitly displayed, but their presence can be directly linked with the six cycle of the full 
matrix map at E — 0. We give a prescription to unravel the other six cycle energy values as 
well, though, E = is the only energy for which we have been able to detect a six cycle of 
the matrix map in the present model. However, in a mixed model of the Fibonacci chain we 
show that following the same prescription we can evaluate at least two values of the electron 
energy at which the matrices will repeat after every six generations. Our method is perfectly 
general, and can be applied to any quasiperiodic sequence. In the latter part of this paper 
we extend the study to the quasiperiodic silver mean chain [10] to show the presence of 
similar kinds of dimers and establish their relationship with the four-cycle matrix map at 
E = already reported in the literature [10]. 

In this context it is worth mentioning that the existence of dimers in random and quasiperi- 
odic chains [13-16] has been known to cause unscatterd/ extended eigenstates leading to a 
perfect transmission. However, the Fibonacci chain and all such lattices in the golden mean 
class, generated for example, using the rule A — > A n B, B — > A [10] have so far remained 
un-noticed in this respect. In a very recent work on a more involved model of the Fibonacci 
lattice Satija et al [16] have shown that there exist a single 'resonant' eigenstate that can 
be explained in terms of the presence of a 'dimer' We begin by studying the oldest and 
the simplest version of the transfer model of a Fibonacci chain. It will be shown that the 
dimers in this system exist naturally, but they combine in a nested fashion to give rise to 
a 'resonance' only when certain conditions are satisfied. In what follows we describe the 
method of analysis. 

The transfer model of a Fibonacci chain comprises of two different 'bonds' L and S generated 
according to the prescription L — > LS and S — > L. We attribute three values to the on site 
potentials, viz, e a , ep and , e 7 corresponding to the sites flanked by L — L , L — S, S — L 
bonds respectively. The nearest neighbor hopping integrals are designated by and t$ 
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corresponding to the L and S bonds [9]. Let us now define three types of transfer matrices 
correspondent to the three sites a, f3 and 7, viz 
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The transfer matrix for any I + 1-th generation chain is obtained recursively following the 
equation, 

M l+1 = M^M, (5) 

where, M\ is the matrix for the Z-th generation, obtained by multiplying the above elementary 
matrices M a , Mp and M 7 in proper order. In the above equation, Mi = M a and M 2 = 
M^Mp = M 7/3 (say) [7]. It has been observed in the literature [7] that in the above model 
with e a = ep = e 7 = and tL — 1 and ts = 2, one gets M; +6 = Mi for I > 2. Let us now 
give a simple physical picture of this phenomenon based on the existence of the dimer like 
positional correlations in this model of the Fibonacci chain. 

Consider the second generation chain, i.e. / = 2. This is LS. The product transfer matrix, 
when written explicitly for this segment, is, M 2 = M- M ; .\/„. We have taken the first site 
of this chain to be of type a without the loss of generality. The last site will be of type 7. 
This is obvious, because the segment LS has to be repeated periodically, so that the bond 
appearing after the S will be L. The matrix for / = 8 i.e. M 8 naturally consists of a much 
longer string of M a , M@ and M 7 . Yet, at E = 0, we all know that Mg = M2 and this is 
repeated after every six steps for I > 2. So, to our mind, a natural question should be, how 
do the individual matrices combine so that the bigger and bigger strings of matrices shrink 
back to M 2 at this particular energy ? This leads us to carefully examine the structure of 
the Fibonacci chain. In Fig.l we show a Fibonacci chain in the sixth generation for example. 
We observe that in it (as well as in an infinite chain), the a-sites appear in isolation and 
the (3 — 7 sites appear in pair, as well as in isolation. It is important to observe that every 
isolated (5 — 7 cluster, as well as (5 — 7 pairs side by side, are flanked by two a-sites. This 
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happens everywhere, even in the infinite chain. Now, if the trace of the matrix M 7/ g = M 7 Mg 
becomes zero, then the corresponding energy renders the matrix product M 7/3 M 7/ g identity 
(barring a negative sign) [12-14]. Once the f3 — 7 clusters (joined by the solid curve in Fig.l) 
combine to give us an identity contribution in terms of the matrices, the a sites at the two 
flanks (joined by a dotted curve in Fig.l) automatically form a 'dimer' [12]. If the trace 
of M a now vanishes for the same energy, then the product of two consecutive M a matrices 
also become equal to an identity matrix (once again with a negative sign), so that a local 
cluster of matrices M a M 1 pM 1 pM a becomes equal to an identity matrix. This will happen 
locally throughout the chain. In this way one can see that the entire string of matrices will 
reduce, in 'size', because of the resonance in these two different kinds of dimers (a — a and 
/?7 — /3j) simultaneously. These pairs appear in a nested fashion, and appear at all scales 
of length because of the self-similar character of the lattice. If we get an energy for which 
= M^p = —I, I being the 2X2 identity matrix, then one can verify that, Mi = M/ +6 for 
I > 2. Now, let us see if E — satisfies this criterion. In the transfer model, 

We select e a = ep = e 7 = and ts = 2 and ^ = 1 [7]. E — definitely makes = —J, but 
Tr(M 7/ g) 7^ and hence, 7^ /. The resonance condition discussed above is not satisfied. 
This of course, does not rule out the possibility of having bigger clusters in the formation 
of nested dimers. For example, one can see (Fig. 2a) that there are clusters formed by the 
triplet fl'-fa and ^7 which are distributed in the same manner and can be considered to 
be a 'renormalized' version of the basic clusters formed by ^7 and a respectively. It turns 
out that in the transfer model the clusters responsible for a six cycle of the matrix map 
at E = are the pairs f3 r yaf3 r y- f3 r faf3 r ) flanked at the two extremities by the cluster afl'j. 
That is, '/'/•(. \/- ; .U„.U, ,) = Tr(M l/3 M a ) = at E = 0. Under this condition M, +6 = M, 
for I > 2. Once this is observed the question that arises immediately is, how to identify the 
minimal clusters responsible for a possible six (or other) cycle of the full matrix map for the 
chain. To resolve this issue we get back to the method of real space renormalization scheme, 
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applied to the infinite quasiperiodic chains in the literature [9]. It is important however, to 
appreciate that, the renormalization group method will be used here only to identify the 
minimal clusters responsible for resonance. 

It is known that the on-site potentials and the hopping integrals in the infinite Fibonacci 
chain renormalize to [9]: 

e a (n + 1) = e 7 (n) + [t L {nf + t s (n) 2 ]/(E - e p {n)); e p {n + 1) = e 7 (n) + t s (n) 2 /(E - C/9 (n)); 
e 7 (n + l) 2 = e a (n)+t L (n) 2 /(E-ep(n));t L (n + l) = t L (n)t s (n)/(E - e p (n)) and t s (n+l) = 
t L (n). 

One defines a parameter space of reduced dimensionality, and defined by the dimensionless 
quantities w, x, y and z, where, 

w n = (E - e a (n))/t L (n),x n = (E - e (n))/t L (n),y n = (E - e 7 (n))/£ L (n) and, z n = 
t s {n) /t L {n). These parameters, under renormalization give rise to the following set of re- 
cursion relations [9], 

w n +\ = {x n y n )/z n — z n — 1/ z n ; x n+ i = (x n y n )/z n — z n ;y n+ i = {w n x n )/z n — 1/ z n , and, 

From the above recursion relations one can show that 

w n+3 = w n+2 w n+ i - w n (7) 

which is equivalent to the trace map (3) with w n = 2x n . Using the set of recursion relations 
(7) together with Eq.(6) it is straightforward to show that if we set w n = w n+ i = then (i) 
w n+3 = w n , and (ii) (w n+6 , x n+6 , y n+6 , z n+% ) = (w n , x n , y n , z n ). That is, the entire parameter 
space is mapped onto itself after six steps of renormalization. Additionally, one comes across 
a three cycle of the trace map (as w n+ s = w n ). As the total transfer matrix across any 
arbitrary generation can be completely expressed in terms of the parameters w, x, y and z, 
it is obvious that the transfer matrix for the generation n + 6 will be identical to that in the 
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generation n for the energy value extracted by setting w n = w n+ \ = 0. Most interestingly, it 
can be checked that w n is actually Tr(M Qj „) and w n+ i = Tr(M Q) „ +1 ) = Tr(M 7/3i „) , where 
Mi jU is the transfer matrix for the %— th site in the n— step renormalized lattice. Therefore, 
setting w n = w n+ i = is equivalent to searching for simultaneous zeros for the polynomial 
equations Tr(M a , n ) = and Tr(M 1 p, n ) = 0. If we start with an infinite Fibonacci chain 
and re-scale it, we observe that the resonance condition, as described above, is satisfied 
for the first time on a two times renormalized version of the original Fibonacci chain (i.e. 
n — 2) if we take E — 0. That is, for E = 0, we get w 2 = and u> 3 = simultaneously. 
The a — a and the 7/3 — 7/3 pairs at this length scale form the required 'dimers' causing 
resonance. To see what minimal clusters do they correspond to in terms of the original 
chain, we trace back to the lattice at the original length scale to see that the elementary 
blocks causing resonance are clusters like a(5^ and P r yaP'y respectively. The 'dimers' in 
a one-step renormalized lattice and their correspondence with the 'dimers' in the original, 
un-renormalized lattice is illustrated in Fig. 2. Doing one more step is simple. It should be 
noted that the method outlined here indicates the possible existence of six cycle for other 
values of E. For example, one could search for simultaneous zero's of, say, w 4 and w 5 etc. 
However, we have been unable to detect any other value for the energy in this model that 
may lead to other cyclic behaviors. But, we propose a second model which clearly shows 
that based on the same thread of reasoning one can get more than one energy values for 
which a six cycle map for the transfer matrix exists. 

A mixed model: 

Let us select a model in which e a = \/5, = e 7 = 0, = 1 and t$ = 2. With this choice 
of the parameters we find that, Tr(M a ) = Tr(M lt3 ) = in the original (un-renormalized) 
lattice if we set E = This leads to the fact that if we start with the above mixed model 
Fibonacci chain at the 8-th generation, then the a — a and P r y—P'j pairs in the original lattice 
will form the nested dimers. The corresponding transfer matrices will combine together to 
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give rise to an identity matrix, and one is finally left with a product of M 7j gM Q in the 
8— th generation corresponding to the triplet of sites aj3^ which is the same as the second 
generation chain. That is, we achieve the result M 8 = M 2 , M 9 = M 3 and so on. Let us now 
renormalize the lattice once and look for a simultaneous solution of the equations W2 = 
and w 3 = 0. Most interestingly we find that w 2 = w 3 = for E = —\^E. If we trace 
back to the original lattice it is easy to discern that the fundamental clusters responsible for 
resonance in this case are a(3^ and ^7 respectively. This implies that for E = —\^E the six 
cycle feature Mi = M; +6 sets in for Z > 3. Thus we have been able to identify two different 
energy eigenvalues for this mixed model of the Fibonacci chain for which the matrix map 
has a six cycle. For the two different energy eigenvalues the 'dimer correlations' are revealed 
at two different scales of length. This is an important feature, linked with the inherent self 
similarity of the lattice concerned. 

The silver mean lattice: 

To check the validity of our arguments further, we have extended this idea to the transfer 
model of the silver mean lattice, which is generated according to the rule L — > LLS and 
S — > L. It may be mentioned at this point, that the entire discussion made for the silver 
mean case is valid, in general, for the class of lattices grown according to the rule L — > L m S, 
S — > L. The recursion relations depicting the parameter space in the silver mean case are 

[17], 

w n +\ = (w n x n — l)y n / z n ~ w n z n — x n j z n \ x n+ i = (w n x n — l)y n / z n — w n z n ; y n+ i = (w n x n — 
l)^n/ %n •K n / z n , z n j r \ {w n x n 1)/ z n , and, t n ^.\ x n y n / z n z n \j z n 
From the above relations it follows that [17], w n+ 2 = w n+ it n+ 2 — w n and, t n+ 2 = w n w n+ i — 
t n+ ±. The basic clusters responsible for dimer correlation in this case are a— a and aflj— aflj 
respectively. We find that, (w n+4 , x n+4 , y n+4 , z n+4 ) = (w n , x n , y n , z n ) if we set w n = w n+1 = 0. 
This happens for E — 0. That means the four cycle of the matrix map [10] is observed 
following our prescription. The trace map exhibits a two cycle behavior in this 
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w n+ 2 = w n . One can look for other possible cycles of Mi by renormalizing the lattice and 
following the arguments given earlier. These aspects are now being investigated and the 
results will be published in due course. 
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Figure Captions: 



Fig.l: Sixth generation Fibonacci chain with three different sites. Clusters forming dimers 
are displayed by solid (for (5 — 7) and dashed (for a) lines respectively. 

Fig. 2: Portion of an infinite Fibonacci chain. Clusters responsible for dimer-like correlations 
in one-step renormalized lattice (b) are the (3 — 7 doublet (dashed box) and the single a sites 
(solid box). They correspond to the triplet f3^a (dashed box) and the doublet (3^ (solid 
box) respectivelyin the un- renormalized lattice (a). 
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